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Analyses of collective flow and space-time evolution based on relativistic 
hydrodynamical model 

C. Nonaka a , N. Sasaki a , S. Muroya b and 0. Miyamura a 

a Department of Physics, Hiroshima University, Higashi-hiroshima, 
Hiroshima 739-8526, Japan 

b Tokuyama Women's College, Tokuyama, Yamaguchi 745-8511, Japan 

We numerically solve fully (3+l)-dimensional relativistic hydrodynamical equation with 
the baryon number conservation law. For realistic initial conditions we adopt the results 
from the event generator (URASiMA). Using this model we discuss collective flow. 

1. INTRODUCTION 

The various kinds of flow phenomena such as directed flow, elliptic flow and radial 
flow have been observed in recent experiment at AGS and SPS Anisotropic flow 

phenomena are expected to show a strong dependence on the nuclear equation of state 
0. Here, for analyzing anisotropic flow phenomena with relativistic hydrodynamical 
model in detail, we need full calculation of (3+l)-dimensional relativistic hydrodynamical 
equation. Our purpose is analysis of collective flow based on (3+l)-dimensional relativistic 
hydrodynamical model with realistic initial conditions. 

2. MODEL DESCRIPTION 



Our calculational procedure is explained briefly. To 
begin with, we calculate initial conditions for the rela- 
tivistic hydrodynamical model by using the event gen- 
erator, URASiMA (Ultra-Relativistic A-A collision 
Simulator based on Multiple scattering Algorithm) 
which is characterized by Multi chain model ||. We 
focus on AGS energy region and several incident en- 
ergies such as 5, 10.8, 20 AGeV for Au+Au collisions 
in which impact parameter is fixed to 6 fm ||. We 
characterize time to when the projectile nuclei finishes 
passing through the target nuclei and hydrodynamical 
evolution starts. Energy density and baryon density 
from URASiMA is transformed into temperature and 
chemical potential with the equation of state of the 
ideal hadron gas. Here, so as to solve the relativistic 
hydrodynamical equation easily we smooth the results 
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Figure 1: The solid line stands 
for freeze-out in the present 
calculation. The dashed line 
and the dotted line are refer- 
enced in 101. 
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keeping condition that the difference between original results and smoothed results is less 
than 1 %. We numerically solve fully (3+l)-dimensional relativistic hydrodynamical equa- 
tion coupled with the baryon number conservation law using Lagrangian hydrodynamics 
[0] . The equation of state is needed to solve the hydrodynamical equation and anisotropic 
flow is evaluated in detail by inputting the different equation of state to the hydrodynam- 
ical equation. Here, for the first trial we use the equation of state of the ideal hadron gas 
which contains resonance up to 2 GeV ||. Furthermore, we assume that the hadronization 
process occurs when the temperature and chemical potential of the volume elements cross 
the boundary (solid line in fig.l). The solid line is so designed that the freeze-out tem- 
perature becomes 140 MeV at vanishing chemical potential, based on chemical freeze-out 
model and thermal freeze-out model 0. Finally, we investigate azimuthal distribution of 
proton by using Cooper-Frye Formula [Kj for particle emission from hadronic fluid. 

3. CALCULATED RESULTS 



t = 5 [fm/c] t = 15 [fm/c] t = 25 [fm/c] 




Figure 2. The evolution of temperature distribution and chemical potential distribution 
at z ~ fm in Au+Au 10.8 AGeV collision. (b = 6 fm) 

Figure 2 displays that the evolution of temperature distribution and chemical potential 
distribution at z ~ fm in Au+Au 10.8 AGeV collision. The graphs at t — 5.0 fm/c show 
initial conditions from URASiM A and from initial time hydrodynamical expansion starts. 
We can see that both temperature distribution and chemical potential distribution expand 
in the transverse direction. The temperature distribution decreases with time, on the other 
hand chemical potential distribution increases with time. This is because baryon number 
density localize and the average of mass of baryons is larger than lGeV. The fluctuations 
which we can see in initial chemical potential distribution descend with time. Freeze- 
out hypersurface which is assumed in previous section is illustrated with temperature 
distribution in fig. 2. The life time and size of fluid depends on the hydrodynamical 
expansion and freeze-out conditions. The size of fluid can be estimated roughly to be 
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about 10 fm at t = 15 fm/c from fig. 2. This value is larger than the experimental result, 
6 fm which is derived by analysis of 2-dimensional tt + it + correlations JTTJ . 



t = 5 [fm/c] 



t = 15 [fm/c] 



t = 25 [fm/c] 
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Figure 3. The evolution of flow at y ~ fm in Au+ Au 10.8 AGeV collision, (b = 6.0 fm) 



Table 1 

The evolution of fluid in different incident energies. 



t — to hyper surf ace 


E Uh [AGeV] Max. e [GeV/fm 3 


Max. n-Q [fm 3 


Max. \v\ 


t [fm/c] 


5.0 1.29 


0.79 


0.19 


8.0 


10.8 2.06 


0.94 


0.28 


5.0 


20.0 2.68 


0.97 


0.39 


4.0 


freeze-out hyp ersur face 


^ [AGeV] Max. e [GeV/fm 3 


Max. n B [fm 3 


Max. \v\ 


Max. life time [fm/c] 


5.0 0.10 


0.083 


0.40 


24.0 


10.8 0.14 


0.091 


0.46 


21.3 


20.0 0.16 


0.083 


0.48 


21.5 



Figure 3 shows the evolution of flow at y ~ fm in Au+Au 10.8 AGeV collision. 
From the initial local velocity distribution we can see that the projectile nuclei finishes 
passing though the target nuclei. Though only small transverse flow exists at initial time, 
it increases with time and the growth of transverse flow depends on initial temperature 
distribution and chemical distribution. We analyze the space-time evolution of the fire 
ball in several initial conditions (table 1). On the t = to hypersurface, maximum energy 
density, maximum baryon number density and maximum local velocity increase with 
incident energy respectively. On the freeze-out hypersurface, Max. e and Max. ub are 
diluted, on the other hand Max. \v\ is larger than initial one. Therefore the internal 
energy density turns into kinetic energy of the flow. The Max. life time of fluid is not 
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Figure 4. The dependence of the trans- 
verse momentum distribution on the 
azimuthal angle in Au+Au 10.8 AGeV 
collision, (b = 6.0 fm) 



proportion to incident energy, because life time of fluid is determined not only by initial 
energy density but also by local velocity on the freeze-out hypersurface. Figure 4 shows 
the Px-azimuthal fluctuation of proton in Au+Au 10.8 AGeV collision. Elliptic flow is 
obviously seen in larger P T region. Detailed analysis will be presented elsewhere [O . 



4. SUMMARY 

By using full (3+l)-dimensional hydrodynamical calculation coupled with the baryon 
number conservation law, we investigate the hydrodynamical evolution of hot and dense 
matter by inputting the results of the event generator, URASiMA for the realistic initial 
conditions. From analysis of the evolution of fluid with this model, we obtain several 
results at AGS energy region. As a result of our freeze-out conditions, the life time of 
fluid is comparatively long and the size of fluid which is estimated roughly is larger than 
experimental result which is obtained by analysis based on HBT ||11|| . By comparing 
physical properties on the t = t hypersurface to ones on the freeze-out hypersurface, the 
decrease of energy density reflects on the growth of local velocity. Finally, from fig.4 the 
influence of elliptic flow to the particle distribution increases with transverse momentum. 
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